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Table 1: Parameters for ensembles processed so far; j3 is the gauge coupling, a is the lattice spacing, L and 
T are the spatial and time-extent of the lattice, respectively, m-^ is the pion mass. 


1. Introduction 

This talk presented our efforts towards a charmed meson physics program on RBC/UKQCD’s 
Nf = 2 + \ domain wall ensembles using the same discretisation for the light and the charm quarks. 
Compared to lattice QCD for light quark physics, simulations for charm are comparatively scarce 
(cf. the FLAG report [1]). One of the reasons is that the charm quark sets a mass scale which 
until recently was very difficult to reconcile with the energy scale set by the light quarks in a fully 
relativistic and dynamical lattice setup. In fact, covering both energy scales at the same time still 
poses a challenge due to the large lattice size that one needs to simulate in order to keep both finite 
volume and cutoff effects reliably under control. 

Our programme consists of first studying domain wall fermions as a charm quark discretisa¬ 
tion within the quenched approximation. Let us emphasise that the intention is not to make any 
phenomenologically relevant predictions from within this framework but rather to understand the 
behaviour of domain wall fermions for heavy quarks and the continuum limit in detail. To this end 
we created ensembles and measurements for mesonic matrix elements for cutoffs in the range of 
2-6GeV. While a similar study with dynamical domain wall quarks is still out of question for the 
foreseeable future we expect a number of properties to be very similar. In particular, we expect 
that the continuum limit scaling observed in the quenched theory over a large range of lattice spac- 
ings will be qualitatively the same also beyond the quenched approximation. We can therefore use 
the scaling behaviour found in the quenched theory to constrain contributions to the scaling from 
higher orders in the lattice spacing when analysing data in the dynamical case but on a smaller 
range of lattice cutoffs. Our exploratory quenched studies have been covered in Tsang’s and Cho’s 
talks at this conference [2, 3]. The main result is the observation of a large region of heavy quark 
masses for which the automatic 0(a)-improvement of domain wall quarks is maintained and spec¬ 
tral quantities and matrix elements extrapolate to the continuum limit linearly in a^. We found that 
by tuning the domain wall height, a free parameter in this discretisation, a reduction in the size of 
cutoff effects can be achieved. The experience gained in the quenched case is now being applied to 
a charm phenomenology program using dynamical gauge field ensembles. This talk reports on the 
status of these simulations. 

2. Our strategy 

Table 1 gives an overview of RBC/UKQCD’s domain wall ensembles which we are currently 
analysing for this charm study. The ensembles to the left have a lattice cutoff of L7-L8GeV 
(coarse) and the ones to the right 2.4GeV (medium). We are currently generating an additional en- 
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semble aiming at a third lattice spacing with k, 2.8GeV and with a pion mass around 200MeV 
(fine). The global strategy can be summarised as follows: we will compute heavy-light, heavy- 
strange and heavy-heavy meson observables (masses, leptonic and semileptonic decay matrix ele¬ 
ments, bag parameters) on all ensembles and also the vector two-point function relevant for com¬ 
puting the anomalous magnetic moment of the muon. Based on the findings in fhe quenched sfudy 
we are cerfain fhaf predicfions for fhese quanfifies for fhe physical charm quark mass and even 
heavier, fowards fhe bottom quark mass, can be made on fhe medium and fine ensembles wifh 
small cutoff effecls and funclioning aufomafic C?(a)-improvemenf. This may nof be fhe case on fhe 
coarse ensemble where fhe range in quark masses where we expecf to have a good confrol over cuf- 
off effecls doesn’f allow fo simulafe direcfly af charm. On fhe coarse ensemble fhe bare inpul quark 
mass corresponding lo charm is above 0.45 and hence in fhe pari of parameler space where our 
numerical evidence from fhe quenched Iheory suggesls fhe beakdown of 0(a)-improvemenf [2]. 
We fherefore consider a slrafegy where resulfs for all quanfifies wifh heavy quark masses slighfly 
smaller fhan charm are made in fhe conlinuum limif. This is Ihen followed by a small exfrapolalion 
in fhe heavy quark mass fowards fhe physical mass of charm. This does nof preclude considering 
a global analysis ansafz faking advanfage of fhe resulfs for heavyier quark masses on fhe medium 
and fine ensemble once fhe latter has been generaled and analysed. 

In fhe following we presenf firsf resulfs for fhe heavy-lighl and heavy-slrange decay conslanf. 


3. Choice of parameters 

Our gauge held ensembles represenl QCD wifh Nf = 2+\ dynamical havours af fwo differenl 
lattice spacings [4, 5, 6, 7]. The basic parameters of fhese ensembles are listed in Table 1. For 
fhe lattices wifh physical lighf quark masses we use fhe Iwasaki gauge acfion [8, 9] and fhe do¬ 
main wall fermion acfion [10, 11] wifh fhe Mbbius-kernel [12, 13, 14]. For fhe unphysical quark 
mass ensembles we use fhe standard Shamir-kernel. The difference befween bofh kernels in our 
implemenfafion [7] is a rescaling such fhaf Mbbius domain wall fermions are equivalenf to Shamir 
domain wall fermions af fwice fhe exfension in fhe hffh dimesion. Mbbius domain wall fermions 
are hence cheaper to simulafe while providing fhe same level of lattice chiral symmefry. Resulfs 
from bofh formulafions of domain wall fermions he on fhe same scaling frajecfory fowards fhe 
confinuum limif where cutoff-effecfs are of 0{a^). All lattices entering fhe following analysis have 
values of mj^L > 4. We fherefore expecf hnife volume effecls to be af fhe percenl level. 

We use Mbbius Domain Wall fermions also for fhe discrefisafion of charm quarks where fhe 
choice of simulafion pai'amefers in fhe dynamical sfudy is based on our hndings from fhe quenched 
case. In parficular, fhe choice of fhe Domain Wall Heighf furns ouf fo be crucial in reducing 
discrefisafion effecls. We also found fhaf discrefisafion effecls become sizeable abruplly as fhe 
bare inpul quark mass is chosen larger fhan aiuh ~ 0.45. This is also reflecled in an increase in 
fhe residual mass. The residual mass is a measure for fhe amounl of chiral symmefry violalion of 
domain wall fermions due fo fhe hnife exlenl of fhe 5lh dimension. The axial Ward idenlily for 
domain wall fermion takes the form 

= 2m(P(x)P(0)) +2(75,WP(0)) . (3.1) 
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Figure 1: Motivated by the lattice axial Ward Identity the observable in the top panel is constant towards the 
centre of the lattice for the range of simulation parameters where the domain wall mechanism is functioning 
well. We found that the residual mass starts increasing severely as the bare input quark mass is chosen 
beyond 0.45. From this point on it’s no longer guaranteed that chiral symmetry is maintained. The bottom 
plot shows the mass dependence of the residual mass (we took the value of at T /2). 


Here, is the conserved axial current, P is the pseudo scalar density and J^q is also a pseudo¬ 
scalar current but it mediates between the boundaries of the 5th dimension and the mid-planes. We 
define the residual mass as 


W5,{x)pm 

_ X 

~ I(P(x)P(0)) 


X 


The behaviour of the residual mass as we increase the bare input quark mass is shown in figure 1 . 
For bare quark masses below around anih = 0.45 we see a plafeau which is fhe expecfed behaviour. 
For larger inpuf quark masses fhe behaviour changes drastically indicafing fhe breakdown of fhe 
domain wall mechanism. Based on fhese findings all simulafions shown in fhe following will be 
for a number of inpuf charm quark masses anih ^ 0.45. The heaviesf charm quark we simulafe 
fherefore corresponds fo masses of 2.8GeV on fhe coarse ensemble and 3.3GeV on fhe medium 
ensemble (for comparison fhe physical rjc mass 2.9836(7)GeV [15]). 

On fhe coarser physical poinf ensemble we generafe 48 time-plane complex Z 2 noise sources 
per configuration and on fhe finer ensemble 32 lime-planes. Meson fwopoinf funclions are fhen 
computed using fhe one-end-lrick [16, 17]. On fhe unphysical poinf ensemble fhe number of source 
planes varies and simulations are sfill ongoing. One parficular concern was fhe convergence of fhe 
conjugafe gradienf minimiser in fhe compulation of fhe quark propagator. We employ fhe lime-slice 
residual [18] 


r(t) 


|Dt/r-T]| 


(3.3) 


where D is fhe Dirac operafor, y/ fhe solution vecfor and rj fhe source veclor. The norm | • |f is 
reslricled to fhe lime-slice t. Global residuals will nol notice convergence issues al large Euclidean 
lime-dislances from fhe source in a siluafion where \Y\t decays wilh t over more orders of magni- 
lude lhan covered by fhe numerical precision. The lime-slice residual is a local (in time) residual 
lhal allows lo conlrol convergence more reliably. 
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Figure 2: Left: data for the D-tneson decay constant vs. the rjc mass. The blue and red bands represent 
the result of a polynomial fit to the data. The dashed vertical lines indicate the reference Tj^-masses which 
we interpolated to. Right: The results for the decay constant obtained for the reference r]c masses are then 
interpolated in the pion mass to its physical value. This interpolation is over a very small range given that 
two of our ensembles have near-physical simulation parameters. 




Figure 3; Left: Preliminary result for the continuum extrapolation of the data with physical light quark 
mass at the reference heavy-quark mass points. Right: Extrapolation of the data in the continuum limit 
(black squares) to the physical charm point (solid vertical line to the right). 

4. Preliminary results 

The idea is to compute heavy-light and heavy-strange meson observables for a number of 
unphysical heavy quark masses below and above the physical charm quark mass as illustrated in the 
l.h.s plot in figure 2 for the case of the D-meson decay constant. These will then be interpolated to 
common values for the T]c-mass on each ensemble (the analysis discussed in the following uses the 
reference masses =2.3GeV, 2.5GeV, 2.6GeV and 2.76GeV as indicated by the dashed vertical 
lines). The mass dependence is close to linear and interpolating with polynomials of various orders 
leads to very similar results. Results from different orders for this interpolation can at a later stage 
be used to devise an estimate of a systematic error which will likely turn out to be very small. We 
note that at this stage we do not yet have full statistics and error bars for all ensembles will be much 
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smaller once the simulations are finalised. 

The r.h.s. plot of figure 2 shows fhe resulfs of fhe above inferpolafion ploffed againsf fhe pion 
mass. The poinfs fo fhe very leff in fhis plof are from our physical poinf ensembles. Alfhough fhese 
ensembles represenf QCD wifh valence and sea pions close fo fhe physical poinf we sfill wish fo 
sfudy inferpolafion ansafze fhaf allow us fo use also fhe resulfs obfained on fhe heavier sea quark 
mass ensembles. Firsfly, fhis will allow for correcting small misfunings in fhe lighf quark mass 
towards fhe acfual physical poinf. Secondly, in fhe near ferm we will nof be able fo generafe fine 
ensembles wifh physical sea quarks. Under fhe reasonable assumpfion fhaf fhe lighf quark mass 
dependence will suffer only very small cutoff effecfs our precise knowledge of fhe mass dependence 
on fhe coarse and medium ensembles will allow for exfrapolaling fhe resulf on fhe fine ensemble 
towards physical lighf quark masses. The inferpolafion shown in fhe r.h.s. plof in figure 2 are done 
under fhe assumpfion fhaf fhe slope wifh is independenf of fhe mass cutoff. Wifh heller slalislics 
on fhe coarse and medium ensemble we will soon be able fo make more rigid slalemenls aboul fhe 
validily of fhis assumption. 

The interpolation in fhe lighf quark mass is followed by fhe conlinuum exfrapolalion. The 
l.h.s. plof in figure 3 can only be indicative - a continuum limil wifh only fwo poinfs for charm 
observables is very risky and no quanlifalive phenomenological conclusions should be drawn from 
our dafa al fhe momenl. This will have fo waif unfil we have resulfs for fhe finer lalfice spacing. The 
r.h.s. plof in figure 3 shows fhe exfrapolalion of fhe resulfs in fhe continuum limil fo fhe physical 
charm poinf (solid vertical line) using a linear ansalz. 


5. Conclusions and outlook 

In fhis falk we presenfed fhe sfalus of our new charm physics program based on RBC/UKQCD’s 
Nf = 2 + \ domain wall fermion ensembles wifh an added valence domain domain wall charm 
quark. We have shown firsf resulfs for fhe D-meson decay conslanf. As anlicipaled based on fhe 
dafa from our exploralory quenched sfudy we expecl ralher mild culoff effecls for our choice of 
simulalion dafa and judging from fhe dafa on our coarse and medium lattice fhis remains so also in 
fhe dynamical case. 

We are currenlly generaling an ensemble for a Ihird lattice spacing which will allow for reliably 
exfrapolaling our resulfs fo fhe continuum limif. Togefher wifh fhe dafa for fhe lepfonic D- and Ds 
decay conslanf we are also generating heavy-heavy correlalors and also fhe correlation funclions 
relevanl for determining fhe meson bag parameler and semi-lepfonic meson decays. Furlher fo 
charm phenomenology we are sludying how fo make use of resulfs for heavier lhan charm quark 
masses which we are generaling in particular on fhe medium and fine ensemble. To fhis end we have 
in mind inferpolaling wifh resulfs in fhe sialic limif (HQET) buf also applying fhe ratio method [19]. 
Acknowledgements: The research leading to fhese resulfs has received funding from fhe European 
Research Council under fhe European Union’s Sevenlh Eramework Programme (EP7/2007-2013) / 
ERC Granf agreemenf 279757. PAB acknowledges support from STEC granls ST/E000458/1 and 
ST/J000329/1 and NG acknowledges fhe STEC granf ST/J000434/1. The aulhors gratefully ac¬ 
knowledge compuling time granled Ihrough fhe STEC funded DiRAC facilify (granfs ST/K005790/1, 
ST/K005804/1, ST/K000411/1, ST/H008845/1). 
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